Current-Induced Step Bending Instability on Vicinal Surfaces 
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We model an apparent instability seen in recent experiments on current induced step bunching on 
Si(lll) surfaces using a generalized 2D BCF model, where adatoms have a diffusion bias parallel to 
the step edges and there is an attachment barrier at the step edge. We find a new linear instability 
with novel step patterns. Monte Carlo simulations on a solid-on-solid model are used to study the 
instability beyond the linear regime. 
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Vicinal crystal surfaces can exhibit a number of differ- 
ent morphological instabilities that may be important in 
crystal growth and device fabrication. One of the most 
interesting of these instabilities arises from the biased 
diffusion of adatoms under the influence of an external 
driving force, as illustrated by the step bunching insta- 
bility on Si(lll) surfaces 0] induced by heating with a 
direct electric current. Many different aspects of this so- 
called electromigration problem have been investigated 
when the driving force is normal to the steps: e.g., the 
microscopic origin of the driving force on adatoms, the 
effect of temperature on the instability, and the resulting 
step patterns ||-||] . We show here that when the relative 
orientation of the external driving force and the miscut 
angle of the vicinal surface is varied, there exists a qual- 
itatively different instability that can be explored both 
experimentally and theoretically. 

Because of the intrinsic anisotropy of a stepped sur- 
face, kinetic instabilities can usually be associated with 
one of the two principal axes normal to and along the 
average step orientation (the x and y directions, respec- 
tively), ft is convenient to describe step positions in 
Fourier space. From a linear stability analysis, the evo- 
lution of small perturbations in the positions of the steps 
away from a uniform step train 5x n (y, t) can be written 
as 



5x n (y,t) =J2 5x (Q^) eiqy+in ' i ' +U ' ( ' q '' t ' )t + 

9,0 



C.C. 



(1) 



where q is a wavenumber along the step direction, and 
<fi a phase factor relating consecutive steps. When 
Rew(g, <f>) > 0, a perturbation of that mode will grow 
exponentially. 

One class of instabilities involves the growth of finite q 
perturbations along the y (step) direction, as illustrated 
by the Bales-Zangwill M instability, fn this example it 
is sufficient to take = 0, since this (in-phase) mode is 
maximally unstable and will dominate. Indeed the un- 
derlying (Mullins-Sekerka) instability can be understood 
by considering a single interface. 

Another class of instabilities involves the relative mo- 
tion of neighboring steps, and can usually be understood 



using a one dimensional (fD) model with q = and the 
step pairing mode <fi = n. Examples include the step 
bunchinginstability due to either the Ehrlich-Schwoebel 
barrier |7]]|], impurity adsorption or the usual 

electromigration instability with the current normal to 
the steps. Stoyanov || used a ID model to study the lat- 
ter instability in detail. He showed that a vicinal surface 
can be unstable towards step bunching when adatoms 
have a diffusion bias in the step-down direction and there 
exists an extra barrier for attachment to step edges. The 
step pairing mode (4> — tt) gives the maximum instability 
unless there are very strong repulsions between steps Jf2| 
or significant direct interterrace diffusion [ p^|Jf4[ . 

Quantitative studies of electromigration based on a 
generalized ID step model have shown very good agree- 
ment with many aspects of the experiments [[l3| . How- 
ever, a more general perspective is required to under- 
stand some recent surprising results by Latyshev et 
al. They observed the sudden appearance of step 

bends in certain regions of the surface where the step 
edges are nearly parallel to the electric current, leading 
to the formation of anti-step (down-step) bunches as well 
as the usual (up-) step bunches. Motivated by this re- 
sult, we study here a 2D step model with biased diffusion 
parallel to the step edge direction, and indeed find an in- 
stability. This new instability requires both finite q and 
finite </> and reaches its maximum when the phase shift 
<f> = ±7r/2 at small q, unlike any other surface instability 
we are aware of. 

Figure [I] gives a schematic view of the surface and 
introduces the notation. As in Stoyanov's model ||, 
adatoms have a drift velocity D s F/kT, where D s is the 
adatom diffusion rate on flat terraces. F is the mag- 
nitude of the external force acting on adatoms, and we 
assume in this case it is in the y direction. For simplicity, 
here we ignore evaporation and deposition of adatoms on 
terraces; this does not affect the basic instability we find. 
We make the usual quasistatic approximation jf6), set- 
ting dc/dt = in the diffusion equation describing the 
adatom concentration field c(x, y) on the terraces: 
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FIG. 1. Schematic view of a vicinal surface with an exter- 
nal force parallel to the step edge. 



To carry out a linear stability analysis [19| , we substi- 
tute the solution of the diffusion equation Eq. (0) with 
boundary conditions Eq. (^) into Eq. ([!]). Keeping terms 
linear in Sx, we find for the growth rate uj(q, 4>)'- 



-2D s c° eq a 2 A q [fqd sin + g{q, </>)] 
2A 9 dcosh(A 9 w ) + (A 2 d 2 + l) sinh(A g w ) : 



(8) 



where 



](q, </>) = g x (l — 2 coscj) + cos 2 0) — g y q 2 cosc/> 
+ [g y q 2 + 2g x (l- coscj))] 
x [cosh(A 9 wo) + A q dsinh(AqWo)] , 



(9) 



\dx 2 



dy 2 



D s Fdc 
kT dy 



0. 



(2) 



The boundary conditions are given by linear kinetic ex- 
pressions with a parameter n governing the rate of ex- 
change of atoms between steps and terraces: 



±«[< q (2/) -c(xt(y),y)] = n- J(x±(y),y), 
where the surface adatom flux J is given by 

J fa, V) = -D a (-d/dx, -d/dy + f) c{x, y), 



(3) 



(4) 



and / = F/kT. The +(— ) indicates the ascending (de- 
scending) terrace next to step n as shown in Fig. [l| 
c n q (?/) i s the local equilibrium adatom concentration of 
step n at position y, and n is the unit vector normal to 
the step edge. Through microscopic mass conservation, 
the normal velocity of step n can be written as 



an • 



[Jfan.y) - Jfa»,y)] 



D s a 2 [2c e M-c{x+, y)-c(x-,y)], 



(5) 



where a is the lattice constant in the x-y plane, assuming 
a square lattice. 

To take into account the stabilizing effects of the line 
tension and step repulsions, we write 



< q (y) 



exp 



kT 



(0) 



where the step edge chemical potential fJ, n {y) — the 
change in free energy per atom for adding atoms to the 
step at position y — is given by |p7|,|l8| 



V'iw^-V'iwn^+^Xn/dy 2 +fi C - (7) 



Here we assume there is a nearest-neighbor repulsive 
interaction V(w) between steps separated by a terrace 
width w and V'(w) = dV/dw. j3 is the step edge stiff- 
ness and nc is the equilibrium atom chemical potential 
of the crystal. 



vV + ifq, 



and 



9x 



d 2 V 



k T dw 2 



W — Wq 



kT' 



(10) 



(11) 



In the above equations there is an important length scale 
d = D s /k, the ratio between the adatom diffusion rate 
and the attachment rate, as introduced by Pimpinelli et 
al. J2TJ]. g(q, (/)) is the stabilizing term from the line ten- 
sion and step repulsions; it reduces to the known result 
of Pimpinelli et al. [ p0| in the case when g x = (i.e., 
ignoring step repulsions) and / — > 0. 

The first term in the square bracket of Eq. (||) is the 
major result of this paper. For / > (external force in 
the +y direction), Re a; is positive for — tt < 4> < and 
is maximally unstable when = — 7r/2 for a particular 
wavenumber q that is not too large [for large enough 
g, the relaxation term g(q,(f>) will dominate]. This is 
quite different from the usual step pairing instability 
with maximum instability at <f> = tt if the bias is in the 
step down direction. Figure ^ is an illustration of a step 
train with <j) — ~n/2. For several neighboring steps, 
the hills and troughs are shifted in their positions and 
form regions of high step densities (narrow terraces) that 
extend along certain direction. Changing the sign of / 
(the direction of the external force) makes modes with 
< < 7r unstable. Note that if the attachment rate is 
very fast (so that k — > oo or d — > 0), the mass transport 
is limited by the diffusion rate and the instability disap- 
pears. Thus the extra barrier for adatom attachment to 
step edges is crucial for the existence of the instability. 

We can gain a better understanding of the physical 
origin of this instability from the following qualitative 
picture. Figure [| draws three steps with perturbations 
in the form of sinusoids with a phase shift <j) — —n/2 
(steps are out-of-phase with every other step). We as- 
sume there is a very strong barrier for adatom attach- 
ment to a step edge. The direction of the atomic flux 
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on a terrace driven by an external force parallel to the 
average step edge direction is then approximately the 
average of the directions of the two steps bounding the 
terrace. 

From microscopic mass conservation [Eq. (^)], the ve- 
locity of a segment of a step is given by the difference of 
the normal components of the two surface atomic fluxes 
at the two sides of the step. For segments A, C, and E in 
Fig. the two fluxes cancel and the velocity vanishes at 
those points. However for segments B and C, there is a 
net contribution from the two fluxes that makes the step 
move in the direction shown. Thus the amplitude of the 
perturbation will increase and the step train is unstable. 
Changing the direction of the external force will reverse 
the direction of the motion and stabilize the perturba- 
tion of Fig. ^, but will destabilize a perturbation with 
phase shift <j) = n/2. Thus a vicinal surface is unstable 
when there is driven diffusion parallel to the step edge 
direction. This picture is only correct when there are no 
efficient channels for adatoms to hop directly from one 
terrace to another one. 

It is natural and straightforward to extend the linear 
stability analysis to an external force in a general direc- 
tion. The general solution is too complicated to write 
down here, but in the limit of weak external forces and 
small wavenumber q, we have 



2d + wq 



1- COS0 

— 2/ a; a—— — ■ f y q dsm< 

2a + wq 



where f 3 



(v) 



(12) 

F x (y)/kT is the x(y) component of the 



(reduced) external force. The term proportional to f x 
has the familiar (1 — cos (f) form for the pairing instability 
for an external force normal to the steps. 

These results may help us understand the origin of the 
anti-step bunches seen by Latyshev et al. on Si(lll) 
at elevated temperatures. Initially the current direction 
is perpendicular to the average step direction and step 
bunching occurs in the usual way. Because of sublima- 



tion, nearly uniform arrays of crossing steps form be- 
tween the step bunches [Q . As the system coarsens and 
the distance between bunches gets larger, the inclination 
angle of the crossing steps increases and becomes nearly 
perpendicular to the step bunches, and thus parallel to 
the current direction. It is at this stage that Latyshev 
et al. [ fl5| observe the appearance of step bends in the 
crossing arrays that grow in size and accumulate into 
an anti-step bunch, i.e., a step bunch in the opposite 
direction to the initial orientation of the vicinal surface. 

Since the physical model discussed in this paper also 
give good agreement with the initial step bunching be- 
havior of Si(lll) |T^] when the current is normal to the 
steps, it seems quite likely that the instability discussed 
above is responsible for the appearance of the step bends 
in the crossing arrays. However, a detailed analysis of 
the anti-step bunch formation in the experiments neces- 
sarily involves the nonlinear regime and a treatment of 
interactions with the preexisting step bunches, which is 
beyond the scope of this letter. 

To get some information about the nonlinear behav- 
ior, we carried out kinetic Monte Carlo simulations of a 
solid-on-solid (SOS) model Jl3| that incorporates the es- 
sential physical features of biased diffusion with attach- 
ment barriers at steps. The external force is mimicked by 
an asymmetry in the attempt frequencies for exchange of 
surface atoms between nearest-neighbor sites. The en- 
ergy barrier is simply the binding energy of the surface 
atom in its initial site, except that there is an extra bar- 
rier when the movement causes a change in total energy, 
and hence for attachment at steps (see Ref. [jL3) for more 
details). Figure || is a snapshot of a simulation with the 
driving force parallel to the average step edge direction. 
The initial step edges are in the vertical direction and 
adatoms have a downward diffusion bias. The dark re- 
gions are step bunches formed by step segments bending 
in the opposite direction to the individual steps (black 
lines) on the terrace. Note that only a relatively small re- 
gion of each step lies within the step bunch, in contrast 
to the usual step bunches that form when the current 
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FIG. 4. Snapshot of a simulation using a solid-on-solid 
model with a diffusion bias parallel to the step edge direc- 
tion. The initial steps are along the vertical (y) direction. 

is normal to the step direction. This general feature is 
already suggested from the patterns in Fig. || 

The earlier SOS simulations of Dobbs and Krug p2| 
gave very different results. They used Metropolis dy- 
namics, which does not incorporate the step edge bar- 
rier, and the resulting step patterns contained several 
artifacts and did not compare well with experiments on 
Si(lll). 

Krug and Dobbs [f2l|]2^ 1 also introduced a continuum 
model that showed that a surface is quite generally un- 
stable under an external driving force when the adatom 
mobility has a slope dependence. This prediction, un- 
like other features, e.g., the orientation of selected facets 
and the coarsening rate, is independent of the underly- 
ing physical origin of the slope dependence (in our case, 
the step edge barrier) and is thus robust. Indeed one 
can make connection between the step model and the 
continuum model by taking <p — > q x wo and q — > q y and 
expanding Eq. ( |l2] ) for small q x (y) : 

LO ~ f x q 2 x + (2d + w Q )w 1 f y q x q y + <D(q A ) (13) 

Except for an anisotropic ratio that is absent from their 
model, this equation is qualitatively the same as the re- 
sult of the continuum model in Ref. |2^] for a surface 
tilted away from a high symmetry plane with the broken 
symmetry of the underlying crystal structure along the 
x direction. Although their continuum model does not 
capture all the essential physics of a vicinal surface below 
the roughening temperature, e.g., the anisotropy in the 
adatom mobility and the singularity in surface stiffness 
at zero slope, it does give an idea of the instability at 
very large length scales. 

In conclusion we report a new linear instability on vici- 
nal surface that arises from the interplay between a step 
train and the two dimensional diffusion fields between 
the steps. There are many interesting theoretical issues 
yet to be resolved. The study of the instability beyond 
the linear regime is needed to understand the coarsening 



behavior and anti-step bunch formation. The effect of a 
current along the step direction is also being investigated 
experimentally p3|. 
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